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Abstract 


Based  on  the  generalized  Levinson’s  theorem,  some  general  properties  of 

quantum  sticking  coefficients  of  an  atom  at  0  K  are  derived.  Good  agreement  is 

4  4 

obtained  with  experimental  results  for  a  He  atom  colliding  with  a  liquid  He 
surface. 


2 


The  problem  of  quantum  sticking  coefficients  a  of  an  atom  colliding  with  a 

solid  surface  has  long  been  a  matter  of  controversy.  Lennard- Jones  and 

Devonshire*  first  calculated  a  based  on  the  distorted-wave  Born  approximation 

(DWBA),  giving  a  vanishing  coefficient  at  T  *  0  K,  or  equivalently,  when  the  wave 

number  k  of  the  incident  atom  goes  to  zero.  The  result,  however,  disagreed  with 

2 

the  results  of  classical  numerical  calculations  and  experiments  both  predicting 
unity  instead  of  zero  for  sticking  coefficients  at  T  +  0.  As  a  possible 

3 

explanation  for  such  a  disagreement  between  DWBA  and  experiments,  Goodman  and 

4  5 

Garcia  and  Ibanez  pointed  out  that  a  Morse  potential  (exponentially  decaying  at 
large  distance)  used  by  Lennard- Jones  and  Devonshire  for  describing  an  adatom- 
surface  interaction  is  unrealistic.  They  considered  an  algebraically  decaying 
potential  and  obtained  finite  sticking  coefficients.  Their  result,  however,  is 
inconclusive  because  the  WKB  approximation  was  erroneously  used  for  the  scattering 
states  of  the  adatom.  Other  explanations  have  looked  at  correlated  motions  of  the 
surface  and  the  adatom.  Krowles  and  Suhl^  presented  a  polaron  picture  that  when 
the  adatom  is  near  the  surface,  it  becomes  heavier  due  to  virtual-phonon  creations 
inside  the  matter.  Although  some  enhancement  was  found  for  sticking 
coefficients,  it  still  turned  out  that  o  -*■  0  as  k  -*•  0.  On  the  other  hand,  Doyen^ 
considered  self-trapping  of  an  atom  near  a  surface  with  a  static  distortion, 
obtaining  the  trivial  results  of  0  <  a  <  1.  Although  the  self-trapping  picture  is 
in  general  more  appropriate  than  the  polaron  picture  for  a  short-range  potential. 
Doyen's  argument  has  a  deficiency:  Doyen  assumed  the  existence  of  a  resonant 
localized  state  for  an  atom  near  the  surface  as  a  quantum  mechanical  analog  of  the 
fact  that  in  classical  theory,  the  atom  strikes  the  surface  with  probability  of 
unity,  but  the  existence  or  nonexistence  of  such  a  state  in  quantum  theory  is  the 
very  core  of  the  issue. 


V- 
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In  this  paper,  we  develop  a  new  theory  of  quantum  sticking  coefficients  at  T 

=*  0  K.  We  have  recently  generalized  Levinson's  theorem  to  an  atom  dynamically 

a 

interacting  with  a  one-dimensional  (1-d)  matter  system.  Based  on  the  theorem,  we 

derive  here  the  following  properties  of  a  at  T  -  0  K:  the  small-k  dependence  of 

o(k)  is  determined  by  the  long-distance  behavior  of  adatom- surf ace  potentials. 

-2 

For  potentials  decaying  faster  than  x  ,  where  x  is  the  adatom- surface  separatior 
a(k)  *  k  for  small  k.  For  potentials  decaying  as  x  or  slower,  a(0)  can  be 
finite.  The  theory  encourages  experiments  on  a  at  small  k 

,-14  4 

<  0.1  A  .  For  a  He  atom  colliding  with  a  liquid  He  surface,  the  present  theory 
predicts  a  critical  wave  number  -  0.01  A  1  below  which  a(k)  depends  on  k 

q 

approximately  linearly,  which  is  in  good  agreement  with  experiment. 

Let  us  consider  the  adatom  +  1-d  matter  Hamiltonian 

Htot  *  H(*»P>  +  V{i’x)  +  K(P)  '  (1) 

where  H(x,p)  is  a  1-d  matter  Hamiltonian  with  x  and  p,  respectively,  representing 
position  and  momentum  vectors  for  atoms  in  the  1-d  matter,  V(x,x)  is  an 
interaction  between  the  adatom  and  1-d  matter,  and  K(p)  is  the  kinetic  energy  of 
adatom.  A  scattering  eigenstate  of  (1)  with  a  wave  number  k  and  a  boundary 
condition  that  !l)(rc,k)  =  0,  where  r  is  an  appropriate  constant  vector  independent 
of  k,  can  be  written  as 

*(r,k)  -  ^  [f(-k)F(r,k)  -  f(k)F(r,-k)]  .  (2) 

Here  r  s  (x,x),  f(k)  =  F(rc,k)  is  a  Jost  function,  and  F(r,k)  is  a  Jost  solution 
having  an  asymptotic  form  e  $Q(x)  41  *  +  ",  where  ♦q(x)  is  the  ground  state  of 
H(JE,p)  (T  -  0  K). 


We  have  recently  shown  that  Levinson's  theorem  in  static  potential 

g 

scattering  applies  to  our  case  in  a  similar  manner: 


Nb  -  6(+0)/ir  , 


where  is  the  number  of  bound  states  of  Ht  ,  and  6(k)  is  a  scattering  phase 
shift  of  the  adatom. 

Based  on  the  generalized  Levinson's  theorem  (3),  we  now  derive  some  general 
properties  of  quantum  sticking  coefficients  a(k),  at  T  =  0  K.  From  time -dependent 
scattering  theory, ^  one  can  write  a(k)  as 


2itM  ^  _ 

a(k)  =  |  6(VEf}  !<flvflk  >1  *  U) 

lc  ^ 

where  ^  means  a  summation  over  final  states  Jf>,  Vf  represents  a  final  channel 

f  + 

interaction,  e.g.,  creation  of  phonons,  and  |k  >  is  the  scattering  eigenstate  (2) 
with  a  trivial  change  of  normalization.  Let  us  consider  k  «  1  such  that  the  S- 
matrix  element  S(k)  =  e2i6(k)  _  ^  and  a  ^arge  distance  x  *  R  such  that 


-U(R)  =  -V(i0,R)/Oi2/2Ma)  -  k2  , 


where  xQ  represents  the  mean  position  of  the  matter  when  the  adatom-matter 
distance  is  large.  Let  us  assume,  for  the  moment,  an  inequality 


kR  <  1  . 


Under  these  conditions  we  have  for  R  <  x  <  1/k 
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»>. 


|k+>  s  -2ik00(x)  .  (7) 

We  now  claim  that  once  (7)  is  realized,  then  for  0  <  x  <  R 

|k+>  «  k  ,  (8) 

where  x  *  0  denotes  the  surface.  This  can  be  seen  as  follows:  over  the  region  R' 

<  x  <  R,  where  R'  is  such  that  k^  <<  -U(R')  «  (-U(x))  ,  |k+>  smoothly  changes 

max 

from  its  behavior  at  x  -  R1  to  its  behavior  at  x  ~  R.  The  behavior  of  |k+>  in 
this  region  must  be  somewhere  between  the  following  two  extreme  cases:  (a)  U(x) 
*  U(R)  and  (b)  U(x)  =  U(R').  But  one  can  easily  see  that  in  both  cases  the 
requirement  of  smoothness  of  |k+>  at  x  =  R  leads  to  |k+>  *  k  at  R'  <  x  <  R.  This 
fact,  in  turn,  leads  to  |k+>  «  k  at  0  <  x  <  R1 ,  because  in  this  region  |k+>  can 
depend  on  k  only  through  a  constant  prefactor  which,  however,  again  by  smoothness 
of  |k+>  at  x  *  R',  must  be  proportional  to  k.  In  this  way  we  reach  (8). 

From  (A)  and  (8)  we  have 

a(k)  «  k  .  (9) 

We  remember,  however,  that  (9)  is  derived  under  the  condition  (6)  and  S(k)  -  1. 

Little  is  known  about  the  small-k  dependence  of  S(k).  Our  analysis  and  numerical 

calculations  below  for  a  Morse  potential  (see  Fig.  1)  show  that  the  latter 

condition  imposes  a  stronger  restriction  on  k  than  (6).  We  now  analyze  (6)  for 

typical  long-distance  behaviors  of  the  potential  U(x).  Let  us  introduce  a 

standard  number  8  <  1  in  analogy  to  the  "standard  value"  in  the  Lindemann  melting 
12 

formula.  With  this  8  we  can  predict  for  each  given  potential  a  critical  wave 


j.v:  .•-V-V-.-J.-.V-V,'. 


number  below  which  a(k)  depends  on  k  linearly  (approximately  when  the  condition 
S(k)  -  1  is  not  fulfilled).  We  have  obtained  the  following  results: 

(i)  For  U(x)  -  -  ye  <x,  a(k)  «  k  at  k  <  k  , 

—  c 

where  -(k  /ic)ln(k^/Y)  ■  p  .  (10a) 

c  c 

(ii)  For  U(x)  -  -  yx  n  (n  >  2),  a(k)  «  k  at  k  <  k  , 

—  c 

where  y^°  k*  =  p  .  (10b) 

c 

(iii)  For  (Ux)  -  -  yx  n  (n  <  2) ,  kc  =  0,  which  strongly 

indicates  a  possibility  of  finite  a(0).  (10c) 


_2 

Numerical  calculations  below  for  a  Morse  potential  (10a)  with  y  “  48  A  and 

<  *  1.4  A  *  (numbers  for  the  W-He  system)  show  that  k£  -  0.1  A  which  leads  to  p 
4  4 

-  0.6.  The  He  atom  colliding  on  a  liquid  He  surface  is  described  by  a  van  der 

Waals  potential  (10b)  with  n  ■  3  and  y  “  20  A.  Using  p  -  0.6  we  have  k£  -  0.01 

"1  9 

A  ,  in  good  agreement  with  experiment. 

We  have  perfomed  a  numerical  calculation  on  a  simple  self -trapping  model. 

The  model  assumes  a  harmonic  lattice,  H  .  =  )  Jta  a+a  ,  a  surface  atom  plus  an 

pn  L  q  q  q 

adatom  interacting  through  a  Morse  potential,  and  the  surface  atom  seeing  another 
Morse  potential  created  by  the  end  atom  of  the  1-d  lattice  in  its  equilibrium 
(taken  to  be  x  ■  0), 


H 


s+a 


P«/2M_  +  0(x  )  +  P^/2M  +  U(x  -x  ) 
s  s  s  a  a  as 


(ID 


The  dynamical  interaction  between  the  harmonic  lattice  and  the  surface  atom  (one- 
phonon  process)  takes  the  form 


In  (12),  is  the  number  of  lattice  atoms,  and  the  unit  step  function  0(t)  means 

that  |k+>  here  is  a  scattering  eigenstate  of  Hg+a.  In  analogy  to  the  second-order 

optical  process,  this  is  a  "luminescence  component"  approximation  thereby 

neglecting  "Raman  components",  which  is  expected  to  be  good  for  treating  short- 

range  potentials  of  the  Morse  type.  We  note  that  Vr  in  (4)  is  H.  at  t  >  0. 

t  mt 

Thus  our  main  task  is  to  solve  the  Schrodinger  equation  Hg+  <J/  *  Eip . 

As  an  example,  we  have  considered  W-W  and  W-He  Morse  potentials: 


,  -2 KX  _  -kxn 
(e  -  2e  ) 


% 


(13) 


where  x^  ~  3  A,  xa  ~  3.6  A,  15  -  1  eV,  D  ~  10  meV  and  k  *  1.4  A  The  one- 

body  problem  for  Morse  potentials  is  exactly  solvable:5  the  W-He  potential  U(x) 
supports  4  bound  states,  whereas  the  W-W  potential  0(x)  supports  about  200,  and 
the  energy  separations  of  the  two  adjacent  lowest  bound  states  are  about  10  meV. 
With  these  facts  in  mind,  we  consider  following  variational  wave  functions  for  the 
eigenstates  of  Hg+a: 


<Kxs,xa)  »  C$q(xs)  (e  a  -  Se  a)  #(xa) 


2  N-l 


*1  l 

n«0  a*l 


C 

pq 


$  (x  )  h  (x  ) 
p  s  q  a 
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(14) 


where  C  =  0  for  bound  states  and  C  =  1  for  |k+>.  Three  $  's  are  confluent 

P 

hypergeometric  functions  representing  the  lowest  3  bound  states  of  the  surface  W 
atom."*  $(x  )  is  an  error  function  preventing  the  adatom  from  penetrating  into  the 
solid.  Finally,  the  h^'s  are  Gaussian  weighted  Hermite  polynomials,^  which  are 
found  to  be  best  suited  for  describing  fully  correlated  motions  of  the  He  atom 
near  the  surface. 

The  Schrodinger  equation  is  thus  converted  to  a  set  of  coupled  algebraic 
equations  for  the  coefficients  S  and  C^'s.  The  stability  and  convergence  of  the 
numerical  calculation  has  been  checked  by  changing  a  certain  set  of  variational 
parameters  as  well  as  the  number  N  of  h^'s.  Figure  1  contains  our  final  results 
for  N  ■  21.  Part  a  shows  |S|  and  ReS  as  functions  of  the  wave  number  k, 
demonstrating  the  generalized  Levinson’s  theorem  (Eq.  (3)).  Part  b  shows  a(k)  vs. 
k.  The  approximate  linear  dependence  of  a  on  k  is  seen  below  k£  »  0.1  A  \ 

To  summarize,  we  have  generalized  the  Levinson's  theorem  to  an  atom 
dynamically  interacting  with  a  1-d  matter,  and  based  on  this  we  have  derived  some 
general  properties  of  quantum  sticking  coefficients  at  T  *  0  K.  A  generality  of 
the  present  theory  strongly  encourages  further  accurate  experiments  on  sticking 
coefficients  in  a  variety  of  physical  systems.  Theoretically,  on  the  other  hand, 
we  must  determine  if  a(Q)  is  finite  in  the  case  U(x)  -  x  ^  at  large  x.  The  result 
of  this  study  along  with  details  of  the  present  paper  will  be  reported  elsewhere. 
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Figure  Caption 

Fig.  1  (a)  | S |  and  ReS  vs.  K.  It  is  seen  that  S  +  1  as  k  +  0.  (b)  a(k)  vs.  k. 
Below  k  -  0.1  A  a  has  an  approximate  linear  dependence  on  k. 
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